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At the mean-field level, on fully connected lattices, several disordered spin models have been shown 
to belong to the universality class of "structural glasses", with a "random first-order transition" 
(RFOT) characterized by a discontinuous jump of the order parameter and no latent heat. However, 
their behavior in finite dimensions is often drastically different, displaying either no glassiness at 
all or a conventional spin-glass transition. We clarify the physical reasons for this phenomenon 
and stress the unusual fragility of the RFOT to short-range fluctuations, associated e.g. with the 
mere existence of a finite number of neighbors. Accordingly, the solution of fully connected models 
is only predictive in very high dimension whereas, despite being also mean-field in character, the 
Bethe approximation provides valuable information on the behavior of finite-dimensional systems. 
We suggest that before embarking on a full-blown account of fluctuations on all scales through 
computer simulation or renormalization-group approach, models for structural glasses should first 
be tested for the effect of short-range fluctuations and we discuss ways to do it. Our results indicate 
that disordered spin models do not appear to pass the test and are therefore questionable models 
for investigating the glass transition in three dimensions. This also highlights how nontrivial is the 
first step of deriving an effective theory for the RFOT phenomenology from a rigorous integration 
over the short-range fluctuations. 



PACS numbers: 

I. INTRODUCTION 

The random first-order transition (RFOT) theorji^ 
of the glass transition builds on a mean-field scenario, 
in which a complex free-energy landscape with an ex- 
ponentially large number of metastable states emerges 
below a critical temperature Td and the configurational 
entropy associated with these metastable states vanishes 
at a lower temperature T^. At Tk a RFOT, i.e. a tran- 
sition with a discontinuous order parameter yet no la- 
tent heat, to an ideal glass takes place. This scenario 
is realized in mean-field-like approximations of a variety 
of glass-forming systems (liquid models, lattice glasses, 
uniformly frustrated systems) as well as in mean-field, 
fully connected, spin models with quenched disorderi^i^ 
{e.g., p-spin model, Potts glass). The latter correspond 
to spin glasses without spin inversion symmetry, and 
their behavior differs from that of the fully connected 
Sherrington-Kirkpatrick Ising spin glass which is charac- 
terized by a continuous transition in place of the RFOT. 
They have been investigated in great detail and have pro- 
vided most of the clues about the generic behavior of 
"mean-field structural glasses" . 

To make progress toward a theory of the glass tran- 
sition, one must however go beyond the mean-field de- 



scription and include the effect of fiuctuations^ in finite- 
dimensional systems with finite-range interactions. The 
main assumption behind the RFOT theory of glass for- 
mation is that the mean-field scenario with a dynamic 
and a static critical temperature retains some validity 
when fiuctuations are taken into account. The ergodicity 
breaking transition at Td is expected to be smeared and 
the metastable states no longer have an infinite lifetime 
because of entropically driven nucleation events. This 
underlies the picture of a glass-forming liquid as a "mo- 
saic state" with its relaxation to equilibrium dominated 
by thermally activated rare events involving "entropic 
droplets" ^ Yet, the main ingredients associated with 
the RFOT are assumed to persist, even if renormalized 
by the effect of the fluctuations^. The issue can be un- 
derstood in the much simpler setting of the liquid-gas 
transition of simple fluids. The mean-fleld van der Waals 
approach predicts a liquid-gas transition with a termi- 
nal critical point. It is known that this homogeneous 
mean-field picture needs to be modified, e.g. through 
the classical nucleation theory and the renormalization 
group: concepts such as metastability and spinodal are 
no longer crisply defined, critical exponents as well as 
nonuniversal quantities are modified, yet the transition 
with the two, liquid and gas, free-energy states remains 
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valid. . . at least when the dimension d is larger than 1. In 
d — 1 the mean-field treatment is plain wrong and pre- 
dicts a transition that is not present. Fluctuations can 
therefore have a more or less dramatic influence on the 
mean-field scenario: this is the key-point concerning the 
relevance of the RFOT theory to glass-forming liquids. 
One natural path to follow in order to investigate the 
effect of the fiuctuations on the RFOT and the two- 
temperature picture is then to consider, both numeri- 
cally and analytically, the various proposed models of 
structural glasses in finite dimensions. The disordered 
spin models are especially convenient as they are both 
well-defined at the mean-field level and much easier to 
investigate than more realistic models of structural glass- 
formers or effective Ginzburg-Landau theories in the 
replica formalism. In particular, they quite directly lend 
themselves to computer simulations and to real-space 
renormalization group (RG) treatments. The major ob- 
stacle on this seemingly straightforward route is that so 
far no traces of the RFOT scenario have been found in 
such studies on finite-dimensional, finite-range models, 
with either a complete absence of transition to a glass 
phase- or a behavior more compatible with a continu- 
ous spin-glass transition^ than a discontinuous, "random 
first-order" , one. 

In this work we clarify the reasons for the discrepancy 
between the solution of mean-field fully connected disor- 
dered spin models and their finite-dimensional behavior. 
All fully connected models that we have studied (Potts 
glasses, M — p spin models and generalizations) provide 
predictions that turn out to be correct in very high di- 
mensions only: because of their infinite connectivity, they 
indeed neglect local fiuctuations and strongly enhance 
frustration compared to their finite-dimensional counter- 
part. 

This stresses the importance of local, short-range fluc- 
tuations on the existence of the RFOT. If one envisages 
integrating out the effect of the fluctuations in a renor- 
malization group setting, it is often possible to proceed 
in a stepwise manner. First, one includes the short-range 
fluctuations up to a given scale and derives in this way an 
effective theory that describes the physics on longer dis- 
tances (or lower energies) . The next stage is then to solve 
the effective theory by accounting for the long-range fluc- 
tuations, a usually highly nontrivial task. In the study of 
critical phenomena and phase transitions, the first step 
is in general bypassed and one relies instead on an effec- 
tive Ginzburg-Landau model based on general symmetry 
arguments or on the mean-field solution of the micro- 
scopic model to which gradient terms are added in order 
to allow for spatial fluctuations. Generically, the type of 
mcan-fleld theory considered does not matter. However, 
it does matter in all of the disordered spin models that 
we have studied. As we shall show, the solution on fully 
connected lattices or in infinite dimension is not repre- 
sentative of the situation in finite-dimensional systems 
because it misses the effect of the (purely) short-range 
fluctuations associated with the mere existence of a lim- 



ited number of neighbors. On the other hand, a better 
account of the latter is provided by the Bethe approxima- 
tion, which, despite being mean-field in character, prop- 
erly describes the local effect of a finite connectivity. The 
case of the Potts glass, which we shall discuss in detail, is 
paradigmatic. In their computer simulation study, Bran- 
gian et al^ found no transition and no glassy behavior for 
a 10-state Potts glass model in three dimensions, contrary 
to what was obtained in infinite dimensions. We show 
here that this is actually expected on physical grounds 
and confirm this through the solution of the model within 
the Bethe approximation. 

The outcome of our study is two- fold: (1) stressing the 
importance to test the infiuence of short-range fluctua- 
tions on models of structural glasses through the Bethe 
approximation or similar approaches and (2) showing 
that disordered spin models mostly appear to fail the 
test and are therefore questionable starting points for in- 
vestigating the glass transition in three dimension. 



II. FROM INFINITE TO FINITE DIMENSION: 
DIFFICULTIES AND PUZZLES 

On the way to finding a RFOT in finite-dimensional 
disordered spin models, several difficulties have been 
identified: 

(1) The RFOT can give way to a continuous transi- 
tion to a spin-glass, with a quite different phenomenology. 
This phenomenon was observed for several models when 
studied by numerical simulations in finite dimensions, see 
e.g. Ref. @]. The lack of evidence for the entropy cri- 
sis predicted by the RFOT theory in finite-dimensional 
disordered spin models has been rationalized by invoking 
the fact that, already at the mean-field level, the relative 
temperature interval between Td and Tk is very small for 
"reasonable values" of the number p of irreducibly inter- 
acting spins in the p-spin glass model or of the number 
q of states in the Potts glass model (say, p,q ^ 10): un- 
der such conditions, the phenomenology should then look 
more like that of a standard Ising spin glass than that 
coming with an RFOT. Eastwood and Wolynes^ have re- 
lated this to the small value of the reduced surface tension 
between glassy states in such models and argued that the 
p's and g's that would be necessary to mimic structural 
glasses to be p ~ 20 and q ~ 1000. 

(2) The RFOT to an ideal glass can be superseded by 
a conventional (usually ferromagnetic) ordering transi- 
tion. This phenomenon takes place for instance in the 
disordered Potts models where a tendency to a simple 
nonglassy ferromagnetic order was observed. Actually, 
this is also what happens in more realistic models of su- 
percooled liquids in which there is always a transition to 
the crystal. In the latter case, however, since the transi- 
tion is first-order one can supercool the liquid and study 
the metastable glassy phase. This in turn is impossible if 
the glass transition is superseded by a second-order phase 
transition, as the above mentioned ferromagnetic one. 
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Either as a result of the above difEcuhies or in try- 
ing to avoid them, it is often found that no remnants 
of the RFOT are present in finite dimensions. Actually, 
it may also be that no transitions at all are observed 
if one has succeeded in getting around points (1) and 
(2). An example is provided by the model studied by 
Brangian et alA. These authors considered a disordered 
Potts glass with a distribution of couplings that is dis- 
placed toward the antiferromagnetic ones and that, as 
a result, has only a small fraction of ferromagnetic cou- 
plings. This is a general procedure to avoid ferromagnetic 
ordering in models without spin inversion symmetry. (An 
alternative suggestion is to study an all-ferromagnetic or 
all-antiferromagnetic Potts model with random permuta- 
tions of the q states; the latter enforce a statistical gauge 
symmetry that prevents ferromagnetic or antiferromag- 
netic ordering^"—.) In practice, Brangian et al. focused 
on a 10-state Potts glass model with a bimodal distribu- 
tion of the couplings having only a fraction x ~ 15% of 
ferromagnetic (positive) ones. Their main result, which 
is often cited as a major problem for the RFOT theory, 
is that the three-dimensional model is not glassy at all, 
despite the fact that the mean-field solution in infinite 
dimension predicts a strong RFOT transition for such a 
large number of states. This, and further analytical and 
numerical work, led for instance Moore et alii to put for- 
ward the drastic proposal that the RFOT never survives 
in three dimensions, its phenomenology being instead re- 
placed by that of a conventional spin glass in a magnetic 
field (which has no transition in three dimensions accord- 
ing to the droplet theoryi^ii^) . 

In the following, we investigate the robustness of the 
RFOT in spin models with quenched disorder that dis- 
play such a transition in their mean-field fully connected 
limit. We clarify how the above listed difficulties arise 
when lowering the dimension from infinity down to three. 
We restrict ourselves to systems with short-range inter- 
actions involving only pairs of nearest-neighbor sites on 
a lattice. This choice is motivated by the fact that such 
models are amenable to real space RG analyses and that 
they are easier to study numerically in comparison to 
systems where multi-site interactions are present, such 
as the p-spin model. We first focus on Potts glasses and 
show that by trying to avoid problems (1) and (2) one 
actually imposes contradictory requirements in finite di- 
mensions: when the distribution of couplings is mostly 
antiferromagnetic, spontaneous ferromagnetic ordering is 
indeed thwarted, but considering large values of the num- 
ber of states q then suppresses frustration and leads to no 
glass transition at all. We then analyze two alternative 
types of models: the M — p spin models introduced in 
Ref. p4] and a new class introduced by us. We discuss 
the physical reasons that make the predictions based on 
the solution of the infinite-dimensional model inadequate 
in low dimensions, pointing out, as already mentioned, 
the importance of short-range fluctuations. 



III. DISORDERED POTTS MODELS 

This section is devoted to a systematic study of dis- 
ordered Potts models on lattices with finite connectivity. 
Our aim is to understand the physical reasons that make 
the discontinuous glass transition (RFOT) so elusive for 
these systems. We shall show that the two confiicting 
requirements discussed above actually put severe con- 
straints on the dimensionality and the connectivity of 
the possible lattices. 

A. Models 

We define in the following the two Potts models that 
we shall focus on: 

(1) Potts glass with q states. Its Hamiltonian is given 

by 

H ^ Jtj 5^,^. , (1) 

where the coupling constants Jij are quenched random 
variables, and the Ui's are N Potts variables that can take 
q different values, {1, 2, , g}. The symbol (i, j) means that 
each pair of nearest- neighbor sites i, j on the lattice is in- 
cluded in the sum only once and ^o-u' is the Kronecker 
symbol. We consider bimodal and Gaussian distributions 
of the coupling constants. The mean-field solution of this 
model was worked out in Ref. and depends only on 
the mean and the variance of the coupling constant dis- 
tribution, denoted Jq and (A J)^ respectively. By scaling 
Jo and (AJ)^ with the spatial dimension d as 

where for convenience we use the rescaled A J as the unit 
of energy and temperature, one obtains a well defined 
mean-field model in the limit c? — oo. (The same is 
true in the limit TV — oo when considering a completely 
connected lattice and replacing d with the number of sites 
N.) 

A detailed presentation of this model is postponed to 
the following section. Here we just recall the main result 
obtained within mean-field theory. Note that we only 
focus on the case g > 4. This corresponds to models dis- 
playing a dynamical ergodicity-breaking transition at a 
temperature Td and a RFOT at a temperature Tk- The 
former is in the class of the singularity found in the mode- 
coupling theory of liquids and the latter is akin to an 
entropy-vanishing, Kauzmann-like, discontinuous glass 
transition. Lower values of q lead to a continuous spin- 
glass transition®. As discussed before, it is important to 
take Jq negative enough in order to avoid ferromagnetic 
ordering. In practice, this means Jq < — ^^2T- 
dynamical and static transition temperatures depend on 
I^Ta^ and Tk ~ for q » li^. The 
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values of the jump in the (overlap or spin-glass) order 
parameter at and both tend to one in the large q 
limit, indicating that the larger the number of states the 
more discontinuous the transition at Tk- 

(2) Random- permutation Potts magnet. Its Hamilto- 
nian reads 

H = J '5o-,,7r,,(a,)> (3) 

<i,j> 

where J > (antifcrromagnet) or J < (ferromagnet) 
and TTij is a random permutation of the q colors that is 
attached to the edge between i and j^. This model has 
a "gauge invariance" that prevents antiferromagnetic or 
ferromagnetic ordering^^^^: indeed, if one permutes the 
states of the spin at a given site i, one can always find 
random permutations associated with all edges emanat- 
ing from i such that the energy does not change; after 
averaging over the quenched disorder {i.e. random per- 
mutations), the (staggered) magnetization is then zero. 
This model has been studied by numerical simulation^ 
in its ferromagnetic version and by the cavity method^*' 
in its antiferromagnetic version. In the latter case, the 
T = limit describes a form of "coloring problem" , the 
states being then interpreted as colors. 



B. Frustration and lack thereof for g S> 1 

In the search for finite-dimensional models displaying 
glassy phenomenology, the disordered Potts models have 
played a central role. They were actually at the root of 
the RFOT theory of the glass transition in the pioneering 
papers by Kirkpatrick, Thirumalai and Wolynesi. For 
this reason and also because the mean-field studies cited 
above predicted a strong discontinuous glass transition 
for a large number of states (or colors) they seemed to 
be the most promising candidates to find a RFOT in 
three dimensions. It thus came as a surprise that simula- 
tion results^ instead showed a complete absence of glassy 
phase for a g = 10 disordered Potts model. 

Our claim, which we substantiate in the following, 
is that the original intuition based on the infinite- 
dimensional limit (or the study of fully connected lat- 
tices) is misguided. A first piece of explanation comes 
by considering the degree of frustration, which is known 
to be central for the existence of glassy phases. This 
is most clearly understood in the case of the disordered 
antiferromagnetic Potts model, which corresponds to a 
coloring problem at zero temperature. For any finite- 
connectivity lattice, such a model becomes unfrustrated 
for q large enough: in physical terms, if the number of 
states, i.e. colors, is too large for a given coordination 
number then it becomes easy to arrange them in a way 
that two neighboring sites do not have the same color, 
hence lifting frustration. Therefore, one does not expect 
any glassy phase for large g on a finite-connectivity, and 
in particular a finite-dimensional, lattice. 



The recent study of the antiferromagnetic Potts model 
with random color permutations on (both regular and 
Erdos-Renyi) random graphs of finite coordination num- 
ber oi^ supports our claim. It was shown that the model 
for 9 > 4 may have a phenomenology similar to that 
of other mean-field structural glasses, with a RFOT ac- 
companied by a dynamical transition at T^. However, in 
agreement with our previous discussion, the RFOT dis- 
appears at fixed coordination number c when the number 
of colors q is sufficiently large: for instance, the transition 
disappears for c < 8 for q = A and for c < 13 there is no 
phase transition already with q = 5. Asymptotically, the 
transition is absent when 2qlog{q) — log(g) — 2 log(2) > c 
(see footnote in Ref. [13 )• For a given lattice coordina- 
tion number, frustration, hence glassy phenomenology, is 
therefore lost when q is too large, typically larger than 
(c/2)/log(c/2). 

The same trend takes place in Euclidean lattices. We 
indeed prove in appendix |^ that the g-state antifer- 
romagnetic Potts model with random color permuta- 
tions has no phase transition on Euclidean lattices when 
q > 2c. This generalizes a proof by Salas and Sokal^^ that 
applies to the standard random antifcrromagnet Potts 
model. In particular we show that there is a unique 
infinite- volume Gibbs measure with exponential decay of 
the correlations at all temperatures and for any realiza- 
tion of the random permutations. Quite generally then, 
frustration vanishes for a large number of colors or states, 
as intuitively expected. 



C. Bethe approximation of the ten-state 
three-dimensional Potts glass model 

In the previous section we have shown that frustra- 
tion vanishes when the number of states q becomes suf- 
ficiently large at fixed coordination number in random 
antiferromagnetic Potts models on both Euclidean lat- 
tices and random graphs. As a result, any putative glassy 
phase transition is wiped out. We now study the effect 
of adding a fraction of ferromagnetic couplings. This 
should introduce frustration and possibly lead to RFOT 
phenomenology. We shall see, however, that this is not 
the case. 

In the following we focus on the model considered by 
Brangian et al."^ , which we have investigated through the 
Bethe approximation: in practice, we have then studied 
the g = 10 Potts glass on a random graph with the same 
coordination number as the cubic lattice, c = 6. We 
have used the cavity method that allows for an analytic 
solution. In this model, the distribution of the couplings 
is bimodal, 

P(J,,) = x5{J,, - J) + (1 - x)5(Jy + J), (4) 

where J = and the fraction of ferromagnetic cou- 
plings X varies from to 1. The value considered by 
Brangian et al. is a; = (2 - ^2)74 ~ 0.146'?.. The phase 
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FIG. 1: Phase diagram of the ten-state three-dimensional 
Potts glass with bimodal disorder obtained within the Bethe 
approximation. The value of x shown by the arrow corre- 
sponds to that studied by Brangian et al.— . We have only 
drawn the transition lines with the paramagnetic state. There 
is also a transition line between the ferromagnetic and the 
spin-glass state which we do not report. 



diagram that we have obtained vi^ithin the Bethe approx- 
imation is reported in Fig. [T] There is no discontinu- 
ous glass transition (RFOT) for any positive fraction x 
of ferromagnetic couplings. Instead, for relatively small 
values of x we find a continuous spin-glass transition and, 
above some threshold, this spin-glass transition is super- 
seded by a ferromagnetic one. The model is not glassy 
at all for a; = (2 - y2)/4 ~ 0.146, at any temperature. 
Although there is no general proof, one expects that frus- 
tration on Euclidean lattices is comparable or less than 
that found on random graphs. Thus, the absence of a 
RFOT on a random graph for x = {2 - ^2) /4 ~ 0.146 is 
fully compatible with the numerical results of Brangian 
et al.. 

From this analysis, we therefore conclude that adding 
a small fraction of positive couplings is not sufficient to 
trigger the appearance of a discontinuous glass transition 
when the connectivity is finite and the number of states 
large. In addition, this shows that the absence of RFOT 
in the three-dimensional model studied by Brangian et 
al. is not primarily due to some long-distance, possibly 
non-perturbative, fluctuations but is just a consequence 
of the local property of the lattice, namely the fact that 
the connectivity is too small compared to the number 
of colors. This effect is correctly captured by the Bethe 
approximation. It is instead completely missed by the 
mean-field analysis based on a fully connected lattice, 
which therefore appears to be quite misleading to predict 
the behavior of three-dimensional systems. 

An additional question that could provide a valuable 
hint for future studies concerns the smallest value of the 
spatial dimension d (and of the associated number of 




FIG. 2: Phase diagram of the five-state Potts glass with bi- 
modal disorder on a regular random graph of coordination 
number c = 18 as a function of x and T, showing Tsp (green 
line), where the paramagnetic phase becomes unstable to- 
wards the antiferromagnetic order, Tiocai (red line), where 
the RS solution becomes unstable towards RSB, and the dy- 
namical temperature Td (blue dotted line). We do not know 
if the virtual coincidence of the points at which Tioc merges 
with Td on the one hand and with Tsp on the other is a result 
of the numerical uncertainty or has a deeper significance. 



states q) such that Potts disordered models show a dis- 
continuous glass transition (RFOT) within the Bethe ap- 
proximation. We consider as a prototypical example the 
case of the Potts model with the distribution of the cou- 
plings given by Eq. (j4]) with q ~ 5 states and placed on a 
regular random graph of coordination number c = 18 
(which would mimic a hypercubic lattice in 9 dimen- 
sions). 

The phase diagram of the model, which we have ob- 
tained through the cavity method, is shown in Fig. [2] as 
a function of x and T. For a small enough x {x < 0.22), 
where most of the couplings are antiferromagnetic, there 
is enough frustration to produce a RFOT, as previously 
shown for x = 0^. This is demonstrated by the fact that 
the dynamical critical temperature, Td, is higher than the 
temperature Tiocai at which the replica-symmetric (RS) 
solution becomes unstable towards replica-symmetry- 
breaking (RSB). However, the glass transition is super- 
seded by an antiferromagnetic transition, which takes 
place at a much higher temperature. We have computed 
the temperature T^p at which the paramagnetic solution 
becomes unstable towards antiferromagnetic order—: see 
Fig. [51 This transition will likely be very difficult to avoid 
in any numerical simulation on Euclidean lattices. As 
the value of x is increased, Tsp decreases, but, at the 
same time, the distance between Tiocai and T^ reduces 
and the glass transition becomes less and less discontinu- 
ous. Finally, when x > 0.22, where there is no instability 
towards the antiferromagnetic order, the glass transition 
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becomes continuous and a conventional spin-glass phase 
is found. 

This example confirms that even in a dimension as high 
as d = 9, avoiding the two problems listed in Sec. |lT] im- 
poses conflicting constraints that are very hard to fulfill. 
Some improvement could come from considering random 
permutations in order to avoid the appearance of an an- 
tiferromagnetic phase. Even in this case, one would have 
to consider for q — i a,t least c > 12, i.e. d > 6 for a 
hypercubic lattice. 



IV. THE APPROACH TO THE 
INFINITE-DIMENSIONAL OR FULLY 
CONNECTED LIMIT 

In this section we study the 1/d expansion for the q- 
state Potts glass on Euclidean hypercubic lattices of co- 
ordination number c = 2d. The aim of this analysis is to 
better understand the regime of validity of the mean-field 
results which, as found in the previous section, appears 
to be rather limited. 

We have derived the 1/d expansion of the replicated 
Gibbs free energy by following the method developed by 
Georges, Mezard and Yedidi a^^'^^ . We have focused on a 
Gaussian distribution of couplings. In order to perform 
the 1/d expansion it is useful to introduce the simplex 
representation of the model, as follows: 



H = - 



q-l 



<i,j> a—1 



(5) 



In this representation q is the number of colors and the 
degrees of freedom Si^a are vectors pointing toward the q 
vertices of a tetrahedron in a {q— l)-dimensional space^^. 
The order parameter within the replica treatment reads: 



(6) 



where a and f3 are replica indices. We refer to Ref. [23| 
for more details on the replica theory for Potts glasses. 
When using the simplex representation, the Hamiltonian 
of the Potts model is very similar to that of the Edwards- 
Anderson (Ising) model. The 1/d expansion can then 
be performed through the high-temperature expansion of 
the replicated Gibbs free energy^!. In the case of Potts 
variables the computation is more cumbersome. For this 
reason we only consider the first order in 1/d. The de- 
tailed computation is presented in appendix B and in the 
following we only report the main results. 

The Gibbs free-energy A can be obtained as an expan- 
sion in Q"^ . We focus on the first three terms only, since 
this is enough to discuss the existence and the properties 
of the glass transition: 



N 



W2 , 

6" 



(7) 



The coefficients of the expansion computed at the first 
order in 1 /d read 



t = l-/3^ 



(g2-10g+10)-2(g-2)/? 



3 "^0 



-2 



Wi = 1 



3^^ 
'2d 



(9-1) 



1112 



(8) 



where 



/3 Jo = /3 Jo + 



(9) 



As in the mean- field solution {d — > oo), we consider a 
one-step replica symmetry broken (1-RSB) form for Q"*^: 
the replicas are grouped in n/m groups of m elements 
and Q"^ is equal to Q for all pairs belonging to the 
same group and to zero otherwise. In order to study 
the glass transition one has to focus on the to -> 1 limit 
of {(5A)/{m — 1) which, within the 1-RSB ansatz, reads 



lim 



>-l TO — 1 



W2 - Wi 3 
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Within this formalism a discontinuous glass transition 
[i.e. a RFOT) is signaled by the sudden change of 
the global minimum of (/3A)/(to — 1) from Q = to 
Q = Qea > 0. A negative sign of the coefficient of the 
cubic term, i.e. W2 > wi, then insures that the transi- 
tion, if it exists, must be discontinuous. In addition, if the 
coefficient of the quadratic term, i/4, becomes negative 
below a temperature Tiocai, then a discontinuous tran- 
sition surely takes place for T > Tiocai- In the strictly 
infinite-dimensional limit, one finds W2 > wi for g > 4 
and Tiocai = 1, in agreement with the exact solution of 
the model^^. (One also has to consider Jq equal to zero or 
negative in order to avoid ferromagnetic ordering.) The 
magnitude of the coefficient of the cubic term increases 
{i.e. the term becomes more negative) with the number 
of states. This fact and the exact mean-field solution 
showing that Qea — > 1 for g — >• oo led to the belief 
that the larger the number of states (or colors) the more 
strongly discontinuous the transition in finite dimension. 
The 1/d corrections also seem to support this conclusion. 
Indeed, at T ~ Tiocai, when d ^ q^, which corresponds 
to the perturbative regime of the 1/d expansion, the cor- 
rections make W2 — wi larger and Tiocai higher when q 
increases (if Jo < 0). The latter condition implies that 
the scaled mean of the couplings, Jo, is negative enough, 
i.e. that the model is antiferromagnetic enough. This 
is however quite surprising since for a Jq that is very 
negative the corresponding finite-dimensional model has 
almost all of the couplings antiferromagnetic. As we have 
shown in the previous section, increasing the number of 
Potts states leads to a decrease of frustration for such 
a system, and eventually a disappearance of RFOT, in 
finite dimensions. The 1/d expansion appears instead to 
predict just the opposite! 
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The disagreement between the two approaches (the 
study on random and Euchdean lattices for a large num- 
ber of states and the 1/d expansion) can be traced back 
to the fact that in the former case one considers the limit 
of large number of states at fixed dimension whereas in 
the latter the number of dimension tends to infinity at 
fixed, although large, number of states. These two dif- 
ferent ways of taking the limits of large dimension and of 
large number of states do not seem to commute. For in- 
stance, in the limit of large dimensions, the couplings are 
such that the average value is typically much less than 
the standard deviation: Jo/d<$^ A J = \f\Jd if d is large 
enough, no matter what is the value of Jq. In conse- 
quence, the starting point of the 1/d expansion is such 
that almost half of the couplings are positive whereas the 
other half are negative, a situation quite different from 
the finite-dimensional one we would like to describe, i.e. 
a Potts glass in which almost all of the couplings are 
antiferromagnetic. 

In conclusion, the mean-field infinite-dimensional the- 
ory appears to be rather singular. It may still be 
representative of the physics of finite-dimensional Potts 
glasses, but only for a very large number of dimensions, 
and certainly not for the three-dimensional case. (As 
discussed above, we expect that d has to be larger than 
the square of the number of states in order to be in the 
perturbative regime related to mean- field theory.) As a 
consequence, the free energy in Eq. ([7]) with the coeffi- 
cients given by Eq. ([S]) does not provide a correct effec- 
tive model for the three-dimensional system even when 
the first 1/d corrections are included. 

V. ALTERNATIVE MODELS 

In the following we focus on two classes of disordered 
models, one introduced in Ref. [23| and one tailored 
by us, which might provide a better alternative to Potts 
glasses. Unfortunately, as we shall show, the drawbacks 
of the disordered Potts models discussed above apply, at 
least partially, to these alternative models too. 

A. The A/-p-spin model 

The M-p-spin models were first introduced in Ref. [IJ] 
and later generalized and studied in Refs. [lTll25| . M 

Ising spins s'f\ a = 1,2,- ■■ ,M are present on each 
site z of a hypercubic lattice. In this work we consider 
the p ~ 3 case, whose Hamiltonian is given in terms of 
products of 3 spins chosen from the spins in a pair of 
nearest-neighbor sites: 

M M 
{ij) a<P 7 

^j7M)^(7)^(a)^(«^ . (10) 
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FIG. 3: Plot oiV = {/3A)/{m - 1) as a function of the overlap 
for the M — p disordered model with Af = 3 and p = 3 at 
Td (dashed line) and at Tk (continuous line). As discussed 
in the main text and in Appendix C, the exact solution of 
the model lead to two kinds of overlaps (or order parameters) 
Q'^' and Q'^^ We plot the Gibbs free energy as a function 
of Q'^' along the line Q'-^-' — Q^^^ /3; a similar behavior is 
obtained for other choices. 



The couplings are independently distributed quenched 
Gaussian variables with zero mean and variance AJ^. 
Note that the interactions involve three spins but only 
two sites. 

The expansion of the Gibbs free energy in terms of 
Q"^ , analogous to the one described above, was worked 
out in Ref. ^25,] for general values of M and p and d = oo. 
The unique value of M such that W2 > Wi, i. e. leading 
to a discontinuous glass transition (RFOT), is M = 3 for 
p = 3. This model was later studied by Migdal-Kadanoff 
real-space RG in three dimensions by Yeo and Mooreii. 
They concluded that no glass transition is present in fi- 
nite dimension and that the behavior of the system re- 
sembles more that of a model displaying an avoided con- 
tinuous spin-glass transition. There is still no numerical 
simulation testing this prediction (work is in progress^^). 

In order to understand the disagreement between the 
two approaches (infinite-dimensional mean-field and real- 
space RG in c? = 3), we have solved the model on a 
fully connected lattice with the same choice of scaling for 
the variance as in Ref. 25]: AJ^ = (this is equiva- 
lent to take d — oo and a variance of the coupling scal- 
ing as 1/d). More details are given in appendix C. We 
have found that the dynamical MCT-like transition takes 
place at Td ~ 0.5970 and that the static RFOT takes 
place at Tk — 0.5963. The plot of the Gibbs free energy, 
{/3A)/{m — 1), within the 1-RSB ansatz is shown in Fig. 
[3]for T = Td and T = Tk- The most striking fact is the 
very small value of the free energy for the whole range 
of overlap. Actually, in the exact solution of the mod- 
els two different overlaps emerge: Q''^'* = llii^?)'^ 
and Q(2) = Y.i{SfS^)^/N. We have plotted the free en- 
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ergy along the line Q^"^^ = Q^^'^ /?> but a similar behavior 
is obtained for other choices. As we discuss below, the 
smallness of the free-energy is possibly the cause of the 
lack of robustness of the mean-field results as the finite- 
dimensional fiuctuations not taken into account within 
mean-field theory are then overwhelming compared to 
the structure of {f5A)/{m — 1) found in d = oo. 

We have not tried to develop a 1/c? expansion for this 
model because in finite dimension, due to the lack of 
inversion symmetry, Si — >■ —Si, a nonzero value of the 
overlap Qo (in the 1-RSB scheme) emerges even at high 
temperature: an amorphous, but trivial, disordered mag- 
netization profile is induced by the quenched disorder. As 
a consequence, the 1/d expansion is more involved (one 
should make an expansion in Q — (3o)• 
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B. The F-model 

We now consider another class of models. As the previ- 
ous one it is characterized by spins with nearest-neighbor 
interactions. The partition function of the model is de- 
fined as 

^ = E n ^'^SlSl Sf) exp(~/3i?) 

where the spins can take the values —1,0,-1-1 and fx gives 
different weights to the 3'^ states per site; /x is a kind 
of generalized fugacity. We focus on the case where the 
only allowed states are {1, 1, 1}, {1,-1,-1}, {—1, —1, 1}, 
{—1,1,-1} and {0,0,0}, i.e ^ = for all other states. 
Moreover, we introduce P as the ratio between /i(l, 1, 1) 
and /i(0,0, 0) and we assume that the first four allowed 
states are all characterized by the same value of fi. (We 
have considered other variants which have all led to sim- 
ilar results; thus, we just present the simplest one.) 
The Hamiltonian of the model reads 



where the Jij^s are independently distributed Gaussian 
random variables with zero mean and variance AJ^. A 
welcome characteristic of this model is that it has a kind 
of spin- inversion symmetry, if one fiips, say, all Sl,Sf 
the probability measure is invariant (likewise for S} , Sf 
and Sf,Sf). In consequence the average magnetization 
is zero at high temperature, i.e. Qq is equal to zero. 

As for the previous model, we have first obtained the 
development in powers of Q"^ of the Gibbs free energy for 
the completely connected model, corresponding to d = 
oo, by using the scaling AJ^ = 1/N: 
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FIG. 4: Plot oiV = {l3A)/{m - 1) as a function of the overlap 
for the F-model for P = 0.05 at Ta ~ 0.6608 (dashed line) 
and Tk ^ 0.6607 (continuous line). Note the scale of the 
free-energy. 



where {Sf""' S^'^) = Q"^ independently of a, i.e. the 
symmetry between the spins is not broken. The coeSv- 
cients of the expansion are equal to 

t=l- P''{{Sif) , = {{Siff , W2 = {SiS2S^)^ 

where {{Si)^) and {S1S2S3) are on-site spin-spin corre- 
lation functions that can be easily computed: see ap- 
pendix D. Because of the property of the probability 
measure /i that we have chosen, it is easy to show that 
((5i)2) = {S1S2S3). Then, by using that ((5i)2) < 1 
(the equality only applies to the P — 1 case that we dis- 
regard), we obtain that W2 > Wi. In consequence, this 
model, as the previous ones, undergoes a discontinuous 
glass transition (RFOT) at a temperature higher than 
Tiocai (always defined as the temperature at which the 
coefficient t becomes negative). 

We have also worked out the complete mean-field solu- 
tion of the model and found that if P is too small the glass 
transition is pre-empted by a standard first-order tran- 
sition from the high-temperature paramagnetic (liquid) 
state directly to the amorphous glass state. For values 
of P larger than P ~ 0.035 the model displays a RFOT 
which becomes less discontinuous as P is increased (the 
barrier at the transition decreases). For instance, for 
P = 0.05 we find Td ~ 0.6608 and Tk ^ 0.6607 and 
Q{Tk) 0.073. As found for the M-p-model, the Gibbs 
free energy is very small for the whole range of overlap 
Q between the first and the second minimum: see Fig. 
ID This feature is observed for any choice of P and for 
all variants of the weights that we have analyzed. It sug- 
gests, as for the Af-p- model, that the predictions of the 
infinite-dimensional mean-field theory remain valid for 
very high dimensions only. 

In addition, we have indeed solved the model on a ran- 
dom graph of connectivity c = 6 by using the cavity 
method. As anticipated, we only find only a continuous 
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spin-glass transition as the temperature is lowered. 



C. Surface tension and a possible origin for the 
lack of robustness of infinite-dimensional mean-field 
glassy phenomenology 

In the infinite dimensional versions of the M — p-spin 
and F-models, we have found that the typical scale of 
the Gibbs free energy between the first and the second 
minimum is very small, much smaller than Tk- In the 
following, we argue along lines similar to those in Ref. 
that this makes the mean-field predictions very fragile, 
except in very high dimensions. 

The Gibbs free energy as a function of the overlap is 
the starting point to compute the so-called "surface ten- 
sion" Y] the latter is equal to the reduction, per unit 
area, of the total configurational entropy that is caused 
by fixing the overlap at the value of the secondary min- 
imum outside a ball of radius R. This quantity plays 
a crucial role in establishing that the glass transition in 
finite dimension is associated to a growing length-scale, 
called "point-to-set correlation length", and a growing 
time-scale^.. Moreover, it is also the starting point for RG 
analyses of the glass transition^i21. Studies performed by 
using instanton technique o^^i^^ have shown that the first 
non-perturbative corrections to mean-field theory lead to 
a value of Y equal to jl^"^ dQ^I3A/{m- 1). Thus, a 
Gibbs free energy that is very small between the two min- 
ima Q — and Q — Qea, as we found for the previous 
models, leads to F ^ Tk- This is problematic because, 
as it has become clear in recent years, some of the fiuc- 
tuations not considered within mean-field theory corre- 
spond to adding an effective quenched disorder and lead 
to local fiuctuations of the surface tension Y and of the 
configurational entropy density S(~2r— . Actually, this is 
expected on intuitive grounds: on small length scales, 
the values of Y and Sc must fiuctuate in an amorphous 
system, especially that envisioned as a "mosaic state" ^. 
The effective quenched disorder affecting Y on micro- 
scopic length scales is expected to be of the order of the 
typical energy scale, i.e. Tk, since it is due to local mi- 
croscopic fiuctuations. In consequence, on small length 
scales the disorder is much larger than the surface ten- 
sion Y when the latter is much smaller than Tk- In this 
case, heuristic arguments, as well as RG-based ones^^, 
suggest that the RFOT found within mean-field theory 
is destroyed. This is similar to what happens for the 
random-field Ising model in finite dimension if the disor- 
der is too strong compared to the ferromagnetic coupling 
(which plays the role of the surface tension Y) . 

The same phenomenon is present in the Potts glass 
model for a reasonable number (less than 10) of states. 
Eastwood and Wolynes^ have argued that one could cure 
the problem and make the surface tension Y large by 
considering a large number of states, q ~ 1000, thereby 
preserving the RFOT in three dimensions. However, this 
is not a proper resolution, as one then encounters an- 



other obstacle: as we have previously discussed, one in- 
deed expects that frustration, and as a result glassiness, 
disappear for large q in three dimensions. 



VI. FRAGILITY OF THE RFOT SCENARIO 
FOR FINITE DIMENSIONAL SYSTEMS AND 
PROPER EFFECTIVE THEORY 

As already stressed, there are two main issues to ad- 
dress in order to understand the effect of fluctuations on 
the RFOT: (1) Once short-range fluctuations are taken 
into account, is the resulting effective theory RFOT-like? 
To be so, the model should favor two phases: one (the 
liquid) in which replica remain uncorrelated and another 
(the ideal glass) in which replica are correlated and dis- 
play a high overlap between typical inter-replica config- 
urations (the latter phase being metastable with respect 
to the former for T > Tk). (2) If the answer to the pre- 
vious question is affirmative, do long-distance, possibly 
non-perturbative, fluctuations destroy the transition and 
alter the scenario developed by mean-field and heuristic 
arguments or not? 

These two issues are mutatis mutandis always present 
in physics but the first one is often of little relevance 
and easily handled in the field of critical phenomena and 
phase transitions. This is not so in the present case. 
We have found in the case of the disordered spin mod- 
els under study a striking fragility of the RFOT sce- 
nario to the effect of short-range fluctuations. Despite 
the fact that the description obtained from the infinite- 
dimensional/fully connected limit displays a ubiquitous 
RFOT, the models instead show in finite dimensions no 
glassy behavior at all or one that is characteristic of a 
conventional spin-glass transition. The answer to the 
above first question is therefore negative for these mod- 
els. This leads us to propose that a model with a putative 
RFOT should better be first tested for the effect of short- 
range fluctuations only, before carrying out a full-blown 
calculation including fluctuations on all scales either by 
computer simulation or RG treatment. This makes sense 
for a schematic model that is itself a caricature of the 
glass-forming liquids one is aiming at describing. For a 
realistic glass-forming liquid model, the question is rather 
whether the RFOT approach has any chance to describe 
its glassy behavior. The issue is then whether one can 
arrive, by some reasonably well-defined coarse-graining 
or RG procedure that integrates the short-range fluctua- 
tions, at an effective theory in the expected universality 
class of the RFOT. We discuss below these two aspects 
without dwelling too much on the first one which we have 
already addressed. 

(1) Does a given microscopic schematic model lead to 
RFOT-like behavior once short range fluctuations are in- 
cluded? 

Most schematic models showing a RFOT in the infinite 
dimensional mean-field limit that have been proposed are 
formulated in terms of discrete variables on a lattice. We 
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suggest to use the Bethe approximation as a first test of 
the influence of introducing a finite connectivity in the 
problem. The Bethe approximation is equivalent to solv- 
ing the model on random regular graphs. These have the 
same finite connectivity as the Euclidean lattices that 
they mimic, but of course not the same behavior at long 
distance due to the absence of loops on all scales and the 
essentially tree-like structure. Thus, within the Bethe 
approximation, one is able to take into account at least 
some short-range fiuctuations (one can always improve 
the results by cluster methods) and as a result get a rea- 
sonable description of the local physics. If short-range 
fluctuations completely change the behavior of a model, 
one expects the Bethe approximation to be able to cap- 
ture this effect. 

As an example, we mention again the 10-state Potts 
glass studied by Brangian et alM-. The fact that no glass 
transition was found in the computer simulation study 
on the three-dimensional cubic lattice could in retrospect 
have been predicted from the result of the Bethe approx- 
imation which already shows the disappearance of any 
RFOT for the same connectivity. The approximation is 
known to overlook the effect of long-distance fiuctuations 
and is mean-field in character. However, it does have 
some merit to study the physics at short distance. 

(2) Is the replicated Ginzhurg- Landau free energy pre- 
dicting a strong RFOT (at the mean-field level) represen- 
tative of realistic glass-forming liquid models ? 

The Ginzburg-Landau functional that has been pos- 
tulated to represent the effective theory of structural 
glassesSS is given in a replica formalism by 
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Tr(g(x)3) 



where Q°'^{x) is a local, but coarse-grained, (matrix) or- 
der parameter describing the overlap among replicas and 
the number m of replicas has to be taken to 1, m — ^ 1+, 
to describe the physics associated with an exponentially 
large number of metastable states (above the RFOT). At 
the mean-field level, the glass transition can be made as 
strongly discontinuous as wanted by increasing W2 — Wi. 
(On the other hand, the transition becomes that of a con- 
ventional spin glass in a magnetic field when W2 < wi.) 

However, besides heuristic arguments about 
coarse-graining'^- and guidance from the infinitc- 
dimensional/fuUy-connected limit [see e.g. Eqs. ([7]) and 
pTjl ], there has been no serious derivation of the above 
functional from a proper renormalization step accounting 
for the contribution of the local, short-range fiuctuations 
in a realistic model of glass-forming liquid. As we have 
precisely seen above how fragile to the latter fiuctuations 
is the RFOT scenario in the case of disordered spin 
models, this step should be a prerequisite to validate 
the RFOT theory as a starting point for describing the 



glass transition in real systems. Long-range fluctuations 
can of course destabilize the RFOT even when starting 
from the replicated Ginzburg-Landau functional with a 
strongly discontinuous RFOT at the mean-field level, 
but matters would be much worse if the RFOT is 
already wiped out by short-range fiuctuations. The 
latter possibility is advocated by Moore and coworkersiii 
who suggest that coarse-graining a glass-forming liquid 
model leads to a Ginzburg-Landau functional akin to 
that in Eq. (jl2l) . but with W2 < wi, hence with no 
RFOT at all. 

A possible means to capture short-range fluctuations 
in liquid models is to consider the system within a cavity 
with amorphous boundary conditions. This method has 
recently been proposed^ and applied^ to measure the 
growth of amorphous order in glassy systems. It amounts 
to studying the thermodynamics of a system constrained 
to have all particles (or spins) outside a cavity in the same 
positions (or configuration) as those of a typical equilib- 
rium configuration. By construction, degrees of freedom 
are then integrated out only up to a length-scale which is 
the size of the cavity. Provided that the size of the cav- 
ity is not too large, for instance by working far enough 
from the putative RFOT so that the point-to-set correla- 
tion length is not much bigger than the size of the atoms 
(which in practice is always the case in computer simu- 
lations), such a numerical experiment may then provide 
the information we are after: if the crossover between a 
high overlap with the reference configuration for a small 
cavity size to a small overlap for a large cavity size be- 
comes sharper and takes place at an increasing length 
scale as one decreases the temperature, this is the sign 
that on the probed length scale, the finite-size system be- 
haves according to the predicted RFOT scenario. It may 
still be hard to go from here to a determination of the 
parameters in the effective Ginzburg-Landau functional 
but it is enough to dismiss the alternative scenario with 
W2 < Wl and a (possibly avoided) continuous spin-glass 
transitionii. In the latter case, the system constrained 
in a cavity would display a completely different behavior. 
In particular, in the first regime, when the point to set 
starts to grow, the overlap would always smoothly de- 
crease as a function of the cavity size, without leading to 
a sharp crossover, and this behavior would be observed 
on a small, never increasing length scale for overlaps of 
the order of one. Only the tail at large cavity size would 
display a decrease to zero, characterized by a range grow- 
ing as the spin-glass correlation length: this however is 
besides our goal of probing only the infiuence of the local 
fiuctuations'^^. The published results of computer simu- 
lations of glass-forming liquids with the cavity method^^ 
point in favor of the RFOT scenario, giving credit to the 
effective RFOT functional as a bona fide starting point. 
Again, we reiterate that passing the test concerning the 
robustness with respect to the short-range fluctuations 
still does not guarantee the fate of the scenario at long 
distance when fluctuations on all scales are taken into 
account. 
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VII. CONCLUSION 

In order to clarify the starting point for studying the 
effect of the fluctuations on the discontinuous glass tran- 
sition (RFOT) found in infinite-dimensional (fully con- 
nected) glass models, we have analyzed disordered spin 
models characterized by interactions that couple spins on 
pairs of nearest-neighbor sites on a lattice. (This restric- 
tion to short-range two-site interactions was motivated 
by the fact that such models are amenable to real-space 
RG analyses and that they are easier to study numeri- 
cally than multi-site interaction systems akin to the p- 
spin model.) 

Proceeding along the way, we have disentangled two 
possible causes for destabilizing an RFOT in finite di- 
mensions. The more obvious, which was anticipated, 
is that long-distance, possibly non-perturbative, fluctua- 
tions could wipe out, or alter, the transition. The other 
one is that already accounting for short-range fluctua- 
tions, as the local effect of a finite number of neigh- 
bors, could drastically change the physics of a finite- 
dimensional system compared to its infinite-dimensional 
counterpart. Our present work is an assessment of the 
latter, which shows the fragility of the RFOT scenario 
to the effect of short-range fiuctuations in the case of 
disordered spin models. This highlights how nontrivial 
is the step of deriving an effective theory for the RFOT 
phenomenology from a rigorous integration of the local, 
short-range fluctuations. 

We have identified and studied several mechanisms 
hampering the existence of a RFOT in finite (not too 
high) dimensional disordered spin models. For Potts 
glasses, the two-fold requirement of having (1) mainly 
antiferromagnetic couplings in order to avoid long-range 
ferromagnetic ordering and (2) a large number of states 
(or colors) in order to have a strongly discontinuous glass 
transition imposes that the lattice dimensionality has to 
be high, typically larger than the number of states. Oth- 
erwise, the model becomes unfrustrated and glassiness is 
wiped out. For instance, we have verified by using the 
Bethe approximation that the 10-states disordered Potts 
model studied by Brangian et al<^ is not glassy in three 
dimensions. The other disordered models that we have 
considered, the M — p spin models^i with M = p = 3, 
and the F-model introduced by us, which both display a 
RFOT in the infinite-dimensional limit, are characterized 
by a very small surface tension y. As a consequence, we 
expect finite-dimensional fluctuations to be overwhelm- 
ing compared to Y and to lead, in a RG sense, to a 
vanishing renormalized surface tension on larger length 
scales and, accordingly, to the absence of a glass tran- 
sition. These results and the l/d expansion performed 
for Potts glasses suggest that the perturbative regime 
in l/d where (infinite-dimensional) mean-field results re- 
main predictive is restricted to very high dimensions only. 

We stress that we have nonetheless not found a funda- 
mental principle or a general mechanism forbidding a pri- 
ori the existence of a RFOT and of the associated glassy 



phenomenology in all possible disordered spin models 
involving only interactions between pairs of nearest- 
neighbor sites in realistic dimensions {e.g., d = 3). How- 
ever, it is a fact — actually a puzzling one — that none of 
these models so far proposed display a RFOT beyond the 
large-dimensional limit. The RFOT theory originated 
from the analysis of disordered spin models. Somewhat 
ironically, structural glass physics turns out to lack ro- 
bustness precisely in these models, whereas it appears in- 
stead to be less fragile at least to the effect of short-range 
fluctuations in liquid models^. Other simple statistical 
mechanical models for structural glasses have also been 
proposed, that usually consider interactions which couple 
degrees of freedom on more than two sites, such as pla- 
quette models"^^ and variants (see also Refs. [39l - l42j for 
recent alternative proposals). Models with multi-site in- 
teractions and no quenched disorder such as lattice-glass 
models do show the anticipated glassy behavior—"— but 
they are prone to simple long-range ordering, in particu- 
lar to crystallization. The addition of quenched disorder 
could then possibly prevent the formation of the crys- 
talline phase without altering too much the glassy one^. 
In any case, we suggest that any new proposal should be 
first tested to assess the robustness of the RFOT to the 
effect of short-range fluctuations, by e.g. studying the 
Bethe approximation or similar procedures. The quest 
for a simple statistical mechanical model showing RFOT 
behavior is still open. 
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Appendix A: Proof of the absence of phase 
transition in the antiferromagnetic Potts model with 
random color permutations on Euclidean lattices 

In this appendix we aim at proving, in the case of Eu- 
clidean lattices, the absence of phase transition in the 
antiferromagnetic g-state Potts model with random color 
permutations when the number of colors/states is large 
enough compared to the lattice coordination number. 
To do so, we generalize the the proof given by Salas 
and Sokal^^ of the absence of phase transition for the 
antiferromagnetic Potts model and apply it to the dis- 
ordered model with random color permutations (note 
that the Salas-Sokal proof also applies to a random- 
bond Potts model, provided all couplings are negative, 
i.e. anti- ferromagnetic). The Salas-Sokal demonstration 
makes use of the Dobrushin-Lanford-Ruelle approach to 
the equilibrium statistical mechanics of infinite volume 
classical lattice systems^ and of Dobrushin's uniqueness 
theorem^. They show that the antiferromagnetic q-color 
Potts model on a lattice of maximum coordination num- 
ber c has no phase transition (at any temperature T > 0) 
when q > 2c. Again a large number of colors/states pre- 
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eludes the existence of a transition. 

The antiferromagnetic Potts model with random per- 
mutations on a d-dimensional Euclidean lattice is defined 
by the following hamiltonian 



(Al) 



where J > 0, < i,j > denotes distinct pairs of nearest 
neighbors on the lattice, the spin variable a can take q 
values {1,2,- •• ,q}, Sa.a' is the Kronecker symbol, and 
TTij is a random permutation of the q colors that is at- 
tached to the (oriented) edge between i and j. This 
model has a "gauge invariance" that prevents antiferro- 
magnetic ordering2r2£: indeed, if one permutes the states 
of the spin at a given site i, one can always find random 
permutations associated with all edges emanating from 
i such that the energy does not change; after averaging 
over the quenched disorder (random permutations), the 
(staggered) magnetization is then zero. 

Following the Salas-Sokal development, we introduce 
the conditional probability distribution for a single spin 
Ui at site i with external conditions given by all other 
spins {(Tfc}' with k ^ i. The new ingredient is that this 
probability distribution is defined for a given realization 
TT of the random permutations on all (oriented) edges of 
the infinite- volume lattice. It is given by the following 
Boltzmann-Gibbs measure, 

PM{<yk}') = Z^({(Tfe}')"'exp[-/3J^(5,^,s,,(,,)], 

(A2) 

where the a priori single-spin distributions over the col- 
ors are left implicit, /3 = l/{kBT), and iXik is equal to Tr^fc 
if the edge is oriented from i io k and to ir'^^ otherwise 
[we have used the fact that 5o.fc,7rfci((Ti) = 5^. Note 
that, as the notations tt is used for the random permuta- 
tions, we have chosen P for the conditional probabilities, 
which is then different from the Salas-Sokal notation^^. 

We also introduce the quantity cj^- that measures the 
strength of the direct dependence between the spins on 
sites i and j: 

C^j = sup{^} {5}^^^^5^^vMjd[P^(cri|{cr}),P^(cri|{?})], 

(A3) 

where d[, ] is the distance between the conditional prob- 
ability measures whose definition is irrelevant here^^. It 
is easy to check that cj^ = when z = j or when i and 
j are not nearest neighbors on the lattice. The interest- 
ing situation is therefore when j is one of the c nearest 
neighbors of i, which will be assumed from now on (we 
consider here a lattice of constant coordination number c 
but the reasoning works as well when c is the maximum 
coordination number). 

Salas and SokaU^ have proven a lemma which implies 
that 



< max 



(A4) 



where the functions and fij are defined as 
/y(c^iki) = exp[-/3J(5^^,5f^^.(^^.)] , 



(A5) 



where iiij is defined from above; is the probabil- 

ity distribution at site i in the presence of all its nearest 
neighbors except j. 



(A6) 



where the notation k/i means that k is one of the c near- 
est neighbors of i and the normalization factor is the 
trace over the q distinct states/colors that the spin at 
can take. Finally, p^ifij] and pf[fij] are short-hand no- 
tations for 



ai=l 



(A8) 



Similar expressions hold for [1 — fij] and [/y]- The 
demonstration of the Salas-Sokal lemma requires that fij 
and fij be in the interval [0, 1], which is verified when all 
couplings are antiferromagnetic as considered here. 



Note that by construction [1 
that 

pn^-hi] _ 



pfik 



'[1-/./ 



f^J]+pf[f^A=l,SO 



(A9) 



-1' 



and similarly for the expression with fij . 

The procedure now involves deriving an upper bound 
for pf [1 - /„] and pj^[l - /y ]. From Eqs. (^1 and (1X9)) 
this provides then an upper bound for cfj . 

An upper bound for p^[l — fij] is derived by noting 
first that from the definition of fij in Eq. (lASp . the only 
nonzero contribution to the sum in Eq. (jA7l) comes from 
states such that ai = nij{aj). There is only one such 
state. The next step is to find a bound for the contribu- 
tion of any of the q states ai to the conditional proba- 
bility distribution pf . Since J > 0, the maximum possi- 



ble weight is when ^ 



0. This occurs 



when fji is different from all 7fi/c((Tfc)'s with k any nearest 
neighbor of i different from j. Since there are c — 1 such 
nearest neighbors, there are at most q— (c— 1) states with 
maximum weight. (The number of states is less when in 
the set {7fifc(crfc)} several colors are repeated.) The con- 
tribution o! f all the other states is strictly smaller. As a 
result the weight of any of the q states ai with the con- 
ditional probability distribution p'^ is always less than 
\/[q — (c — 1)]. Putting together the above results we 
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arrive at the conclusion that [1 — fij] < l/[q— (c— 1) 
From Eq. (|X9|) it then follows that 



(AlO) 



and similarly for the ratio involving fij. 

One then concludes that for all pairs of nearest neigh- 
bors i, J on the lattice 



1 



(All) 



For any given realization of the quenched disorder, the 
so-called Dobrushin constant a'^ = supj X^j^i '^Tj there- 
fore satisfies 



a"" < 



q 



(A12) 



The Dobrushin uniqueness theorem^ states that there is 
a unique infinite- volume Gibbs measure when the condi- 
tion a'^ < 1 is satisfied. Under the same condition, an 
additional theorem implies that the correlations in the 
infinite-volume Gibbs measure decay exponentially with 
distance (see Ref. |19|] and references therein). By using 
Eq. (|A12p . the above condition amounts to 



q > 2c. 



(A13) 



This is true for any realization of the random permu- 
tations on the infinite-volume lattice. This proves the 
absence of phase transition and the exponential decay of 
the correlations, hence the absence of any glassy behav- 
ior, when the number of colors is larger than a threshold 
depending on the lattice coordination number. 



Appendix B: The 1/d expansion for the Potts model 

As discussed in the main text, we use the simplex rep- 
resentation of the Potts glass: 



9-1 



H — ^ Jij ^ Si^aSj^a 
<i,j> a=l 



(Bl) 



In this representation q is the number of colors and the 
degrees of freedom Si^a are vectors pointing toward the q 
vertices of a thetrahedron in a (g— l)-dimensional space. 
These vectors satisfy the following relations: 



9-1 



qSs 



(B2) 



(B3) 



s=l 



(B4) 



/q-l \^ q-1 

E«' ^iq-2)J2<ei+q-l. (B5) 

\a=l / o=l 

Integration of the partition function replicated n times 
over the quenched disorder gives: 



= E n 

<i,j> 



P'^ \ ' \ " aa qP ca cP 

2d ^ 2^'^i,a'^i,b'^i:'^^3,b 



(a,/3) a.b 



+ /^^O E E ^ta^ta 



(B6) 



1. Large-dimension (or large-temperature) 
expansion 

The large-dimension expansion is a kind of high- 
temperature expansion, as discussed by Georges and 
Yedidia^^. Following Refs. [2lll22| we introduce a small 
parameter v in front of the terms that couples different 
sites. We also introduce two sets of Lagrange multipliers 
A"^ and -qf to ensure that the following relations hold 
for each value of the small parameter v. 



and 



As a result, we have 
- m{Qf}, {<}; {Af }, m) = lnTr{5=} exp 

+EE^f (e^^X 
+EE''" (E^"a-('?-i)< 

i a \ a J 

where TV is the number of sites in the system and 

rr _ /^^ CQ qfi qa qP 

^ ^2d^ ^ 2^ i^<^ i^b^l''^ hb 
{a, 13) a,b 

+ 13 Jo E E E ^ta^ta ■ 
(ij) a a 

We have implicitly assumed that by symmetry 
Qlab = ^abQf Va 



(B7) 



(B8) 



vH 



(B9) 
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mf Va . 



and by using Eq. (|B7I) we obtain that ~pAi is equal to 



pose A"'' so that 



To preserve the relations in Eqs. (jB7P and (jBSP we im- 



9a: 



a/3 







, which gives Q"'^ = ^ 

(BIO) 

and rjf so that 

0=0, which gives ™^ = -1^ ^ 5,^, . (Bll) 

* a 

A large-temperature/large-dimension expansion of 
—I3A can be obtained expanding exp(u7J) around -u = 
and putting v = 1 at the end^i. The terms of the ex- 
pansion are expressed in terms of averages ( )q with the 
same weight as in Eq. (IB9p . except that v is set to zero: 

{0{{S-})), = 

Tr{sn I exp f E StaSL ^ il ~ ^)Qf] 

{ \_ i a,l3 \ a / 

+ EE'?? (E'^"--('?-ik 



(B12) 



where 

Zo = 



Tr{s=} exp 



EE^f E^?X-(«-i)Qf 



+EE'?" E^"^-(9-ix 

i a \ a / 

The expansion is of the form 

with, at the zeroth order, 
and at the first order, 

d{-fiA) 



dv 



(B13) 



(B14) 



(B15) 



(B16) 



Then, 



= (^) = 1^ E E E (^?a^i:.^"a5f ^ 

(jj) (q,/J) a, 6 

/5^oEEE(^?a5j:,)o (B17) 



2d 



E E +/3'^o(<z-i) E E« 



<i J> (a, ,3) 



<i,j> a 



To more easily compute the second order we introduce 
the operator 



dX, 



a/3 



^ - (^) + E E ^ I E staSt - (9 - l)Qf 



(",/3) 



EE|^fE^?.-('^-i)<) 



which is such that 

d{0) 
dv 



(B18) 



(B19) 



Useful properties of the operator U are listed in Ref. [2l| 
The term A2 can be written as 



■I3A2 



d^i-^A) 



(B20) 



where Uq is the operator U such that all the averages it 
contains are those with u = in the weight. Focusing on 
the m — case (no ferromagnetic ordering) we get 



(B21) 



2 

(a,l3) a,b 
(7,5) c,d 



(a,l3) a.b 
7 c 

i(/3Jo)'E(g-l)gfgf 



as in the last sum the terms with a ~ /3 give only a 
constant. 



2. Zeroth- and first-order terms 



The first-order term is directly obtained as 



2d 



(B22) 



<i,j> {a,p) 
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On the other hand, the zeroth-order term reads 
-/3Ao = 



hiTr{s„} exp 



i (a,l3) \ a / 



(B23) 



In the following we compute this expression in an expan- 
sion in Xf^ and Q"^ ■ First, we pull out the term not 
involved in the trace, which gives 



« (a,/3) 

inTr^snll-pfE^r'^E^eX 



(B24) 



Next, we expand the argument of the trace to cubic order 
in the Ai's: 



inTr,5=}n-pfE^r^E^?x) 

lnnTr{snfl+ E \"'E^?X + 

i \ {a,fS) a 

1 \ " •,q;^-,7i5V^ Qa qP 07 qS , 
n \ / , '-'i.a'-'i.a'-'i.c^i.c + 



(B25) 



(",/3) 
(7,5) 



^ \^ \'^^ \^ Q'' Q<^ Q 



c 



(a,/3) a,c,e / 

From Eqs. (|B2p and (|B4|) it is easily realized that one 
has to pair the replicas to obtain nonzero contributions 
from the trace. As a result, we find 



(a,fS) 

(a,f}) a,c 
(7,<5) 

i^(A"^)2(,-l)g", 

io',13) 



(B26) 



(B27) 



and 



(«,/3) 
(7,5) 

(«^,0 



i ^ (A"'')3(q - l)(q - 2)g- + iTrA3(g - l)g" , 



where we have used Eq. (|B2[) and 



E E<^?^c =(9-l)(g-2)/g^ (B28) 



,.?J,c \s — 1 



Finally, the expression of the zeroth-order contribution 
is obtained to a 0{X'^) as 



'PAo^ N{q-1) i- J2 ^'"''^Q"'^ +n^^<l (B29) 

\ (",/3) 

4 E (^"')' + ^ E (^"')' + ^TrAM . 

(a,/9) (a,,3) / 



3. Second order 



We start with Eq. (jB21[) . To compute this expression 
as a function of A"^ and Q"^ we again expand the weight 
in ( )q around X"^ = to a 0(A*). The only contributions 
that can be nonzero involve 



qQ^7 _ / qa qI3 q1 
'^3 — \'^i,a'^i,b'^i,c 



(B30) 



with a ^ l3, for 

1. 7 7^ a and 77^/? 

2. 7 = a and 7 7^ /3, 
and 



oa qP q7 qS 



(B31) 



(",/3) 



with a ^ 13 and ^ ^ 5, for 

1. 77^a, (Jt^q;, 77^/3 and 5 ^ P 

2. "f — a and (5 7^ /3 

3. 7 = a and S = (3. 

In the end we need an expression that is cubic in the 
Q"^^s. For each term we therefore make an expansion 
up to the third order in X"^ at most. The detailed com- 
putation is quite cumbersome. Here, we only present the 
results term by term. 

Let us start with the first term, S!^^'' = (^S^S^S^ 
with a ^ l3. 

1. When 7 7^ a and 77^/?, the terms of order 0(1) 
and 0(A) are zero. Hence, one only has terms of 
order O(A^) that become of order 0(A*) when the 
square is taken in the Gibbs free energy and can 
therefore be neglected. 
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2. When 7 = a, the only term of order 0(1) is zero 
and we thus have to compute terms of order 0(A) 



and 0(A ). We find that the hnear term S^y' reads 

q-^X°'^J2^aSbSc (B32) 

s 

whereas the quadratic term iS^^a reads 

Ip-'ix'-'rT. (T.SaS?s^^s:\ (y^s^s^x) 

aiM2 \ S / \ S / 

V^a \ S ) 



Now we consider the terms from S2^^^ = (^S^S^SJS^^ 
with a ^ P and 7 7^ (5. 

1. When 7 7^ a, 7 7^ /?, (5 7^ a, and 5 ^ (i, the terms of 
order 0(1) and 0{X) are zero and there is no need 
to compute the term of order O(A^). 

2. When 7 = a and 6 ^ (3, the term of order 0(1) is 
zero but not the term of order 0(A): 



(B33) 



3. For the quadratic term, there are several contri- 
butions that will be denoted according to the way 
spins are grouped in the sums: 

Sf^fiA, 2, 2) = g-iA"^A"* S:S^S^S2^ (B34) 



^^"jf (3,3,2) = q-2^'^5(^a/3 ^ ^ 



(B35) 



01 \ s 



S2^,f (2,3,3) ^^q-^X^'r Sac X 
T.(T.S^Sisi](j:s',Sisi] (B36) 



ai,a2 \ S 



Note that it is not needed to go to the term of order 
0{X^) because 



( )o = (( )o ^ ^achd + SacSbdf 



(B39) 



= (( )o ^ ^ac6bdf + 25ac6bd{{ )o ^ ^achd) + {^acSbdf ■ 

Since 6ac6bd ( )o is a- constant, only the first term 
matters and computing the term ( )g to the order 
O(A^) is enough to give a result of order 0{X^). 

The linear term reads 

(B40) 

For the quadratic one we find several contributions: 



J2 (Y^aS^s:^sA [Ys^bS^dS^A. 



(B41) 



ai-a2 \ S 



5gf (4,2,2) = l{q-l)5aJbd + E(^"')' J 

(B42) 

Sf^f (3, 3, 2) = ^ A"" A"^ X (B43) 



ai \ S 



4. The Gibbs free energy 
In Eq. (jB21| one has to consider the following terms: 



E E(^"^f^c) 



{a,l3) a, 6 
7 c 



/ 



E\ ' / qa qP 07 q^ 
/ , Wa '^b '^c ^0 

(7,5) c,d 



(B44) 



(B45) 



S?lf (2,2,2,2) =y X'^^X^'SacSbd (B37) T4 - ^ ^ <5.AdQ"^Q^' (^."^f 5,^5^)^ . (B46) 



4. When 7 = a and (5 = /?, all terms are nonzero. The 
term of order 0(1) is given by 



{a,P) a,b 
(7,(5) c,d 

To compute these terms we use the relations obtained in 
Ref. [23i] after introducing where 



qapap rp qa qa qP qP ^ s: 

O4 ^ — -LlSaOjj — Ua, 



(B38) 



Vabc — ^E '^a'^fc'^c^ 



(B47) 
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and 



where 



Fabcd ~ SaSbScSd^ 



(B48) 



and using the Einstein convention of summing over re- 
peated indices, one has 

VabcVade = qFbcde - q^SbcSde (B49) 
VabcVabd^q^{q-^)Scd (B50) 
vl, = q\q~-l){q-2) (B51) 



Faabc = q{q - l)<56c 

Faabb = q{q - 1)^ 



(B52) 
(B53) 



FabcdFabef = 9^ ( '^—^Fcdef + (Jcdf^e/ ) (B54) 



Fl,,^^q\q-l){q^ -iq + i) , (B55) 
from which it also follows that 

Sab^cdVabeVcde = . (B56) 



In Eq. (|B44|) . only the terms in which a = 7 or /3 = 7 
are different from zero. This means that 



s3 = EE(^"^f^c 



(B57) 



After substituting the terms of order 0(A) and O(A^), 
one finds 

S3 = EEf'?"'^"'«'^^^ + VMA"^)'x 



a^l3 a,b 



Faca,a,Vba,a2 + 9"' E ^'"''^"^^abc 



(B58) 



which finally gives 

S3 = (g-i)(g-2)^(A"^)= 



(B59) 



^{q -l){q- 2f J2 i^"^)' + 2(9 - - 2)Tr(A3) 



In Eq. (|B45p the nonzero terms are those for which 7 
or S is equal to a or /3 while the other is different from 
both, or such that 7 and 6 are equal to a and /3. Hence 
we have 



S4 — S4 + S4 



(B60) 



oa \ \ / Qa q/3 Qa q6 

^4 — 2-^ 2-^ Wa ^6 "^c 



(B61) 



and 



S4 = ^EE(^a^^''^"5; 

a=^P a,b 
c^d 

Eq. (IB61[) can be reexpressed as 

^4 = E J2i^"'^--^''d + q-'X''^X'''Fa 

a^p a,b \ 



(B62) 



CK^S C, 



g-2A^*(A"^ + A"^)^i;aca,i;6rfai + 

ai 

^q-\XP')Hac J2 ^ba,a,Vda,a, + SacSbd ^ A^"A"M 
aia2 Vt^'^ f 



(B63) 



which gives 



= - 2{q ~ 1)2 - 2(9 - 1)' E(^''')' 

-4(g- l)2Tr(A^) . 



(B64) 



Eq. (IB62[) can be reexpressed as 



s4 = ^Ee(9"'^"'E----^'^- 

c,d 

+ ^^"'(A"^)' Faca,a,Fbda,a,+ 

+ \[q-l)5a,.Kd (E(A''')' + E(^"'')M + 



(B65) 



ri^P 

which gives 



Si = i(<7-l)(g-2)2^(A"^)2+ (B66) 

a^p 

\{q- l){q 2f Y (A"^)' + (9 - - 2)2Tr(A3) . 
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Finally in Eq. (jB46p . we only need the nonzero terms of 
order 0(1) and 0(A), i.e., 



a,b Vct?4/3 



(B67) 



a=^li Q#/3 ai / 



which can be reexpressed as 



T4 = (9 - 1)2 E Q'^'Q'^'^'' + ^('^ - 1)' E 



Pits 



After putting all these results together in Eq. (|B21[) we 
find that —^f3A2 is equal to 



-2(q-l) + i(g-2)2)E(A" 



-2(g-l)2 + i(g-2)-^)E(^ 



2\ rr„/\3\ 



+ (-4(g - 1) + (g - 2)2) Tr(A3) - 2{q - l)Tr{Q^\) 



+ 



+ ^JoiV(g-l)(g-2) 



E(^"')' + ('?-2)E(^"')'+ 



2 Tr(A^ 



i(/3Jo)'diV(<Z-l)E(Q"'')' 



(B68) 



The final piece of information that is needed is the 
relationship between A"'^ and Q°'^ . This results from 
the equation 



= 



(B69) 



which gives 



ga/3 ^x'^f^ + ^{q- 2)(A"^)2 + E A^^A""^ , (B70) 

V 

A"'3 = Q"/3 _ i(g - 2)(Q"^)2 - ^ Q""Q'"3 (B71) 



For the second-order term we therefore have 



--/3A2^Niq-l) 



^Jo{q-2) + ^{(3Jo)'dj + 



(B73) 



while for the zeroth-ordcr term we obtain 



/3Ao = N{p - 1) 



-iE(Q"')^ 



^(g-2)E(Q"''f + ^Tr(Q3) 



^2 

and for the first-order term, 



(B74) 



/?^i=7V(g-l)^EW"^) 



(B75) 



After collecting all the pieces together, we finally find 
(formally setting v — I) 



-PAc:iN{q - 1) 



E(Q"' 



^_ 1 
T ~ 4' 



(B76) 



^(g^ - lOg + 10) + ^Mq - 2) + i(/3Jo)2d j + 



Appendix C: Mean-field equations for the M-p spins 
model 



As in the main text, we focus on the M = 3 and p = 3 
model. In this case, by following the standard procedure 
described in Ref. [4^, one obtains a mean-field solution 
in terms of the two overlaps Q*^^-' and Q*^^) introduced in 
the text. Below, we only reproduce the final result. Let 
us first introduce the notation 



and 



(A"^)2 = (g"^)2 -[q- 2)(Q"'5)3 - 2Q"^ E 

(B72) 



(CI) 



exp 



7^Q 
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The self-consistent equations for the two overlaps read: 



n 



8exp 



2, QW+3Q(2) 



(C2) 



)(2) _ 



n 



dzf dz^ {E{s}^(P^Q^'\Q^'\{S''})s's^y 



8 exp 



2 /QW + 3Q(2) 



The free energy [divided by (m — 1)] is given by 



(C3) 



F 



W) + >QW + log(8) 



/ 



n ^^.iog.(/3,Q(^',Q'^',{5n) (C4) 



,2 ,'Q(i) + 3Q(2) 



8 exp 



Appendix D: Mean-field equations for the F-model 

We assume a 1-RSB ansatz for the matrix Q"^ and 
let Q denote the intra-state overlap. The corresponding 



mean-field equations can be derived by usual methods 
and read 



Q 



[E.;M'(.S)S'cxp(/3v/^E„ga^ 



E5M'(5)exp(/32S^^(5J2) 



(Dl) 

where 5 is a short-hand notation for Si, S2, S3 and 

The effective 



= exp 



field u has to be determined self-consistently from the 
equation 



2 
2 



1/2 



277 



a=l 



(D2) 



1 E5M(^)(^i)^exp(/3^^Ea(^a)^) 
'2 E5M'5)exp(/32u^J5,)2) 



By solving numerically these equations we have obtained 
the results described in the main text. In the high- 
temperature, paramagnetic, region regime where (5 = 0, 
on-site spin averages can be obtained by just using the 
measure /x'(iS). In this way, one can compute the corre- 
lation functions entering in the expansion of the Gibbs 
free energy. 
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